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From the scattering of semicoherent-state wave packets at high magnetic field, we derive analyt- 
ically the transmission coefficient of electrons in graphene in the quantum Hall regime through a 
smooth constriction described by a quadratic saddle-point electrostatic potential. We find anoma- 
lous half-quantized conductance steps that are rounded by a backscattering amplitude related to 
the curvature of the potential. Furthermore, the conductance in graphene breaks particle-hole sym- 
metry in cases where the saddle-point potential is itself asymmetric in space. Fhese results have 
implications both for the interpretation of split-gate transport experiments, and for the derivation 
of quantum percolation models for graphene. 
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The quantum point-contact geometry formed by 
metallic split gates in the quantum Hall regime is a 
cornerstone of many experiments in two-dimensional 
electron gases (2DEGs) based on semiconducting het- 
erostructures. For example, in recent interferometry ex- 
periments made with electrons, the quantum point con- 
tact plays the role of an electronic beam splitter X At the 
theoretical side, the consideration of a smooth constric- 
tion represents a simple toy model of elaborated quan- 
tum transport theories. For instance, the transmission 
coefficient through a saddle-point potential^ is a central 
piece of the percolation network models,—^ which have 
been introduced to describe the inter-plateaus dissipa- 
tive transport in the quantum Hall effect. 

Smooth constrictions are usually modeled by the local 
potential profile at the bottleneck of the constriction 



V(r) = by 2 - ax 2 



(1) 



where a and b are real positive coefficients characterizing 
the potential (for convenience, we chose the electrostatic 
potential value at the saddle as the origin for the en- 
ergies). Important insights on the tunneling processes 
determining the transport properties through quantum 
point contacts can be gained from the quantum me- 
chanical motion of the electron in such a simple po- 
tential as given by Eq. (|T|). In the standard 2DEGs 
usually described by Schrodingcr's equation, the energy 
dependence of the transmission coefficient through the 
quadratic potential in Eq. (JTJ) is well-known and given 
under high magnetic fields (i.e., neglecting Landau- level 
mixing while keeping the magnetic length Ib finite) by 
the expression^ 



T n (E) = 



1 + CXp — 7T 



E-{n+l/2}(Hu c + () 



(2) 



for the nth Landau level, where Ib = \J ^ c /(\ e \B) is the 
magnetic length, ftw c is the Landau-level spacing in the 
2DEG, and Q = l 2 B {b — a). Tunneling processes give rise 
to a nonzero probability for the electron at energy E < 
to be transmitted on the other side of the constriction 
while an electron at E > goes weakly backscattered 
through the available channels. 

With monolayer graphene consisting of carbon atoms 
packed in a two-dimensional honeycomb lattice appears 
a new host material,— where electrons are confined to 
two dimensions, yet with some exotic properties. The 
observation of an anomalous quantum Hall effect in 
graphene£i£ understood in terms of a relativisticlike spec- 
trum of low-energy electrons^ has been followed by 
numerous experimental and theoretical contributions^ 
aiming at exhibiting specific signatures of the 2D mass- 
less Dirac fermions in a non-uniform potential. The stud- 
ies of simple analytical problems, such as the quantum- 
mechanical motion of massless particles in the quadratic 
saddle-point potential are of valuable interest for 

identifying such unusual properties. However, it turns 
out that quadratic potentials, which are exactly solvable 
in ordinary 2DEGs at any magnetic field, become gener- 
ally not analytically solvable within the Dirac equation. 
Here, we consider the regime of large magnetic fields, 
in which the Landau level spacing in graphene is large 
enough so that one can work in the single Landau-level 
limit. We have recently shownii using a semicoherent- 
state Green's-function formalism that quadratic prob- 
lems then become soluble in this regime. 

The determination of the transmission coefficient 
through a smooth constriction for graphene in the high 
magnetic field limit [i.e., the counterpart of Eq. ([2])] is the 
main result of this Rapid Communication, that we start 
by discussing in relation to conductance quantization in 
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graphcnc, with a detailed derivation using scmicoherent- 
state Green's functions making the rest of the Rapid 
Communication. The exact transmission coefficient for 
graphene in the absence of Landau-level mixing reads in 
the nth Landau level (with n > 1) 
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where e = ±1 is a band index characterizing the electron 
and hole-like contributions, and 



E n ,e = n( + e 



y/ n (hQ c 



C74 



(4) 



with fi c = \f2vpjl-B (vf is the graphene Fermi velocity). 
The transmission probability for the lowest Landau level 
(n = 0) is given by T (E) = 1/2. 

The zero-temperature conductance at chemical poten- 
tial \x is given in terms of the transmission probabilities 
by the Landauer-Buttiker formula 
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where we have accounted for the spin and valley degen- 
eracies in graphene with the overall prefactor 4. 

Then, as found in Ref. [ID for the 2DEG case, the 
conductance quantization in the graphene case shown in 
Fig. [T] directly follows from the transmission probabil- 
ities. The first obvious observation is the half-integer 
quantization of the conductance in terms of the con- 
ductance quantum (here 4e 2 //i) with plateaus at values 
(n + l/2)4e 2 //i, reminiscent of the half-integer quantiza- 
tion of the Hall conductance. Two different configura- 
tions, symmetric and asymmetric with respect to 7r/2 ro- 
tation of the saddle-point potential in Eq. (QJ , have been 
considered, which yield to two different curves. A sym- 
metric saddle-point potential is characterized by a = b, 
thus C = 0, and the conductance shows clear particle- hole 
symmetry with respect to change in the energy sign. In 
the case of an asymmetric saddle-point potential (a ^ b), 
£ becomes non-zero, signaling a breaking of particle-hole 
symmetry in the energy levels Eq. and resulting 

in a non-uniform shift of the conductance steps, which 
is more pronounced for the highest Landau levels. It 
is worth noting that the asymmetry of the potential in 
Eq. ([T]) has a different consequence in the case of the 
standard 2DEG, where it just leads to a redefinition of 
the Landau- level spacing, see Eq. @, and such a small 
quantitative modification appears difficult to perceive in 
an experiment. In graphene, the effect of the potential 
asymmetry should be more easily seen in experiments, 
since it yields an asymmetry between the positive- and 
negative- energy dependences of the conductance. 

Now, we focus on the derivation of Eq. ©. We con- 
sider a single-particle Hamiltonian model for an electron 
of charge e = — |e| and of Fermi velocity vp confined to 
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FIG. 1: (Color online) Zero-temperature conductance for a 
saddle-point electrostatic potential with l%a = l%b = 0.05fi c 
(symmetric case) and l B a = Al 2 B b = 0.1f2 c (asymmetric case). 
Asymmetries of the electrostatic saddle-point potential reflect 
in asymmetries between the electron and hole sectors for the 
conductance. 



a two-dimensional graphene sheet in the plane (x, y) in 
the presence of both a perpendicular uniform magnetic 
field B = Bz and an electrostatic (scalar) potential term 
V{r) [given by Eq. ©] 
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with II = — i/iV r — eA(r)/c, where A(r) is the vector po- 
tential defined by the equation V r x A(r) = B. Here 1 
corresponds to the unity matrix in the pseudo-spin space 
(representing electron and hole degrees of freedom) and 
c is the speed of light. For convenience, we will omit 
both physical spin and valley indices, thus assuming that 
the two valleys of graphene remain completely decoupled 
form each other and can be studied separately. In addi- 
tion, we shall not consider the effects of ripples or a mass 
potential but these could be studied following Ref. [Til 

To describe the electron dynamics at high magnetic 
field, it is usefulii to introduce the graphene vortex states 
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with z = x + iy and Z = X + iY . Here R = (X, Y) is 
a doubly continuous quantum number corresponding to 
the guiding center position in the plane, n is a positive 
integer, A is a band index (defined for a given n), which is 
equal to ±1 if n > 1, and for n = 0. States in Eq. ([7]), 
which can be written as ^„.R.A(r) = (r|n, R, A) within 
the Dirac bracket notation, are eigenstates of Hamilto- 
nian © in absence of the electrostatic potential (V = 0) 
with the energy quantization E n _\ = X^/nhfl c . De- 
spite being nonorthogonal with respect to the degeneracy 
quantum number R 
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the set of quantum numbers |n, R, A) obeys a complete- 
ness relation 
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In fact, the states in Eq. ([7]) form an overcomplete basis 
of states, which have the coherent states character with 
respect to the quantum number R. 

Relation (fTTj) allows one to project the elec- 
tron dynamics onto the vortex representation. We 
can then introduce the vortex Green's function 
G(ni, Ri, Ai, ti; 112, R2, A2, ^2), which gives the probabil- 
ity amplitude for a vortex with circulation (or Landau- 
level index) n\ and band index Ai that is initially at 
position Ri at time t\ to be at point R2 at time £2 with 
a new circulation rii and a band index A2. In the fol- 
lowing, we consider the dynamics projected onto a single 
Landau level, meaning that the vortex circulation is con- 
served (n\ = 7J2 = n). Formally, this corresponds to 
taking the limit vf — > +00. Within a single Landau level 
?i, the retarded Green's function takes the formii 

G n;Ai;A2 (Ri,R 2 ) = (R 1 |R 2 )e^/ 4 ) A — [g n;Xl . X2 (R 13 )] (12) 
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where Ar is the Laplacian operator. In the absence 
of Landau-level mixing and for a quadratic saddle-point 
electrostatic potential with the spatial dependence given 
by Eq. ([1]), the function g can be calculated exactly-^ 
with the result in the energy representation (i.e., after 
Fourier transformation with respect to the time differ- 
ence t\ — t<z = t) 
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where r(t) = (1/ v ItI ) tanh.(y |7|i). Here E is the energy 
and + an infinitesimal positive quantity. The parame- 
ters 7 and C, are geometric coefficients characterizing V(r) 
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The coefficient 7 is directly proportional to the Gaus- 
sian curvature of the electrostatic potential, which turns 
out to be negative for a saddle-shaped (or hyper- 
bolic) quadratic function V(r). Finally, the functions 
^ri;Ai:A 2 (0: which contain the full dependences on the 
Landau-level index n and on the band indices Ai and A2 
read for n > 1 
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and ho ; o;o(t) = e~ ltt >/ 2 for the lowest Landau level n = 0. 
It is worth stressing that the previous expressions are 
valid for Landau-level indices n not too high, for which 

y/n)Ml c holds. 



the inequality |C|, vM < {V n ~+ 1 

The action of the differential operator exp [(/|/4)Ar] 
on the function <?(R) in Eq. (fT2)) is evaluated a*^ 
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Then, inserting Eq. (fT4"]) into Eq. (fT5]l . we can perform 
the Gaussian integrals over u to get 

r+°° _jgit(E+iO+) 
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withf(t) = (l + t(T(t)- 7 T 2 (t)y\ 

We note from Eq. (|16|) that for n > 1 and Q 7^ the 
function h n[ \ i; \ 2 (t) is not diagonal in the A space, indicat- 
ing that A is gcnerically no more a good quantum number 
in the presence of asymmetric saddle-point potentials. A 
straightforward diagonalization shows that e = ± ap- 
pears instead as a good number, with h n - c (t) = e - ltE ^,^ 
According to the above formula (|4]) for E njC , the quantum 
number e clearly labels the electron-like and hole-like en- 
ergy bands. We shall henceforth represent the Green's 
function in this e representation, where the latter takes 
a diagonal form. 

In order to determine the transmission coefficient with 
a given energy channel (i.e., at Landau- level index n 
and band index e fixed), we only need the Green's func- 
tion G n;e (Ri, R2) when the states at vortex positions Ri 
and R2 correspond to the same energy and are asymp- 
totically far from the saddle point located at the ori- 
gin. For a saddle-point potential of the form (TTJ, this 
means taking the limits \Xi\ — > 00 and IX2I — > 00, 
while V(Xy,Y x ) = V(X 2 ,Y 2 ) = const. Making the 
change in variable s = d [l — |7|r 2 (i)] / [l + |7|r 2 (i)], 
where d = \X1Y2 — ^Yil/Zg, we can easily take the 
limit d — > +00 in the integral in Eq. (|19p , and obtain the 
expression 
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where a = sgn^!^) and T(z) is the Gamma func- 
tion. At infinity the vortex is close to the asymptotes 
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In the valence band the vortices carry a positive charge. 
Therefore, the correspondence between the transmission 
(or reflection) process of a negative charge carrier and 
the sign of a is now reversed, so that we have 



i exp 
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FIG. 2: (Color online) Schematic to identify the Green's func- 
tions for reflection and transmission. Semiclassically the tra- 
jectories of vortex states lie on contours of constant potential. 
Solid (dashed) lines correspond to an equipotential line V > 
(V < 0) and therefore for E > to a vortex in the conduction 
(valence) band, respectively. Dotted lines denote the asymp- 
totes of the potential in Eq. fT]). 



of the saddle-point potential V(R), i.e., Y± ~ X\ y/a/b 
and Y 2 ~ —X^yfafb. The two modes e = ± circulate in 
the same direction, and are not mixed since they are well 
separated in energy. 

Within the present saddle-point geometry, the trans- 
mission coefficient can be extracted from the asymptotic 
form of the retarded vortex Green's function, thus per- 
forming scattering theory in terms of coherent-state wave 
packets instead of the more standard plane waves. In the 
conduction band (e = +), the transmission of a vortex 
from the left half-plane to the right half-plane is described 
by Green's function expression [Eq. f|20[) ] with the sign 
function a = — while the vortex reflection, where the 
vortex remains in the left half-plane, is characterized by 
a = + (sec Fig. [2]). Therefore, the ratio of the trans- 
mission amplitude to the reflection amplitude for a given 
Landau level n in the conduction band is given by the 
ratio of these two Green's functions: 
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The relations between the transmission probabilities T n t 
and the transmission and reflection amplitudes, i.e., 
|^n,e| 2 = 1 — l r n.e| 2 = T n e , finally provide the result in 
Eq. ©. 

For the peculiar case of the lowest Landau level, 
contributions both from the original conduction and 
valence bands arise. We get for the electron-like excita- 

tions T+(E) = [l + cxp(-7r(£-C/2)/yH 
and for the hole-like excitations T^(E) 

1 + exp (n{E — C/2)/vTtIJ . Summing up these 
two contributions and considering the equipartition of 
the current between the two types of excitations yielding 
a 1/2 prefactor, we finally get the already quoted lowest 
Landau-level contribution 



T (E) = -[T+(E)+T Q -(E)] 



(23) 



In conclusion, we have calculated the transmission co- 
efficient for a quadratic saddle-point electrostatic po- 
tential in graphene, and found that shape asymmetries 
generic to quantum point contacts break particle-hole 
symmetry in the conductance. Our results should be 
relevant for future split gate experiments in graphene as 
well as for the formulation of quantum network models. 
Theoretically, we have also presented an alternative way 
of deriving transmission coefficients from the scattering 
of coherent-state wave packets. 



1 I. Neder, N. Ofek, Y. Chung, M. Heiblum, D. Mahalu, and 
V. Umansky, Nature 448, 333 (2007). 

2 H. A. Fertig and B. I. Halperin, Phys. Rev. B 36, 7969 

(1987) . 

3 J. T. Chalker and P. D. Coddington, J. Phys. C 21, 2665 

(1988) . 

4 B. Kramer, T. Ohtsuki, and S. Kettemann, Phys. Rep. 
417, 211 (2005). 

5 This result for the 2DEG can be derived from Eqs. (1.2)- 
(1.6) of Ref. in the limit ui c — > oo by neglecting the 
Landau- level mixing correction to the curvature energy E\. 
Such corrections lead to additional yet small contributions 
to the energy dependence of the transmission coefficient in 
the large cyclotron frequency limit, without affecting its 
qualitative behavior. 

6 A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 



(2007). 

7 K. S. Novoselov et al, Nature (London) 438, 197 (2005). 

8 Y. Zhang, Y. W. Tan, H. L. Stormer, and P. Kim, Nature 
(London) 438, 201 (2005). 

9 V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 
146801 (2005). 

A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009) . 

1 T. Champel and S. Florens, Phys. Rev. B 82, 045421 

(2010) . 

2 M. Biittiker, Phys. Rev. B 41, 7906 (1990). 

3 T. Champel and S. Florens, Phys. Rev. B 80, 125322 
(2009); ibid 80, 161311(R) (2009). 



